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^ ■ Abstract 



We discuss several methods of calculating the DIS structure functions F 2 (x, Q 2 ) based 
on BFKL-type small x resummations. Taking into account new HERA data ranging 
down to small x and low Q 2 , the pure leading order BFKL-based approach is excluded. 
Other methods based on high energy factorization are closer to conventional renormal- 
ization group equations. Despite several difficulties and ambiguities in combining the 
renormalization group equations with small x resummed terms, we find that a fit to the 
current data is hardly feasible, since the data in the low Q 2 region are not as steep as 
the BFKL formalism predicts. Thus we conclude that deviations from the (successful) 
renormalization group approach towards summing up logarithms in 1/x are disfavoured 
by experiment. 



1 Introduction 

Recent measurements of the proton's structure functions [JXJ, |2|, |3| have raised the question 
whether the observed rise of F2 at small values of x is due to resummation effects described 
by the BFKL equation [|J. In contrast, this observed rise has already been predicted 
using the conventional renormalization group equation (RGE) a la Altarelli-Parisi f| and 
evolving from a low starting scale || 0. That even data at low Q 2 can be described by 
the RGE is now widely accepted [§, |9| . So the problem remains if one can find a unique 
signature in the current F 2 -data for BFKL- or RGE-based evolution equations, or if there 
is even the possibility of combining the different approaches, aiming at a unified evolution 
equation, which might lead to a better description of the structure function data. 

The first question has already been answered in pi. An update with recent data 



appeared in (T^j. These calculations were based on the approach of Askew et al. fT2|| , 
treating the BFKL-equation as an evolution equation in In 1/x, and calculating the struc- 
ture functions from the resulting unintegrated gluon distribution via the fc^-factorization 
theorem [JTBL ITJ]]. With a consistent choice of the infrared parameters, we obtained steep 
structure functions incompatible with the data. 

In section ^| we have a look at the high energy factorization that was developed by 
Catani et al. 0, [14|, 15| . This corresponds to the small x limit, since s ^> Q 2 means 



x ~ Q 2 /s <C 1. Based on this work, Forshaw et al. [|I(J calculated F2 and Fl by only 
considering the gluon sector, and by only taking into account small x data (x < 1CT 2 ). 
We will show the limitations of this method in section [| 

In section ^ we try to modify the conventional RGE with the complete small x re- 
summed anomalous dimensions of Catani et al., which include also the quark sector in 
next-to-leading order. In a first attempt Ellis et al. |T7| evolved existing parton distri- 



butions with this approach. Blumlein et al. [18] later focused on the implementation of 



energy-momentum conservation. Employing these ideas, we will try to fit the current data 



with special emphasis on the low Q 2 region while constraining the partons in the large 



x region. We supply the details of our previously published results |LT| and extend the 
analysis of the conserving factors. 



Our results confirm earlier related work by Ball et al. [19], POfl , in which energy- 
momentum conservation was implemented differently There three small x parameters 
of the input distributions were fitted to experimental data and the large x parameters 
were matched to existing conventional parton distributions. We note that a recent re- 
summed calculation by Thorne |H[] using leading order "physical anomalous dimensions" 
P^] apparently comes to contrary results. We will summarize our results in section [^, 
and comment on the limitations of BFKL inspired calculations of F% at relatively low 
momentum transfers Q 2 . 

2 High energy kj- factorization 



Catani et al. 14, HJ proposed a method of summing small x logarithms which is based on 



the high /^-factorization. It has the advantage of being compatible with the conventional 
collinear factorization. We just want to present the basic outline, without going into 
details. First, we define the renormalization group equations in Mellin space: 

dfi(n,Q 2 ) ^ 2 

j 
Our interest now focuses on the quark singlet combination f s = xY2i=i(Qi + 9i) an d the 
gluon density f 9 = xg. Note that we have rescaled the parton distribution by x for 
convenience. This implies that the anomalous dimensions are connected to the x-space 
splitting functions Vij via 

7*i(n) = I x n Vij{x) dx = V i5 {n + 1), (2) 

Jo 

which differs by one from the standard definition. These anomalous dimensions are in 
conventional perturbation theory given as a power series in a s , of which the first two 
orders are known today. 



With the method of Catani et al., we are now looking for corrections to the anomalous 
dimensions and coefficient functions of the form a l s +k /n k . Up to now, these have been 
calculated for i — to all orders in a s /n, and partially for i = 1. The singularities in 
moment space (n — > 0) correspond to the small x limit in rr-space. 

The corrections to the gluon anomalous dimension can be inferred from the solution 
of the characteristic equation (a s = Caocs/^', for QCD Ca = 3, Cf = 4/3, Tr = 1/2): 

a s x(lL)=n; x{l) = 2^(1) - ^(t) ~ ^i 1 ~ l)- ( 3 ) 

The gluonic small x resummed anomalous dimensions can then be expressed as 

Igg = 1L\ l gq = yflL- (4) 

The quark sector only contributes in the next-to-leading order, introducing an ad- 
ditional factorization scheme dependence in the expressions. A simple way to find this 
contribution is to look at the physical observable F 2 . One can calculate its scaling viola- 
tions in n-moment space in the limit n — > from the gluon distribution via the following 
expression ((e 2 ) denotes the mean squared electric charge of the involved quarks): 



dR 



?_ _ l<?\(n.9^ X9AU \ f9 
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()hx(f r){Cfr„ + 2N / >y„)F (5) 

= (e 2 )h 2 ( lL )R n ( lL )P (6) 

The function hi is basically the Mellin-transform of the off-shell cross section da/dlnQ 2 , 
whereas R n is a process-independent, but (in general) factorization scheme dependent 
renormalization factor. They are given by [|14| 



W 7 ) = ^1 N T ^ 2(2 + 3 7 -3 7 2 ) / 7T 2 7 2 1 (7) 

2n 3 — 2 7 \1 — 4 7 2 sin(7T7) tan(7T7) ' 



and 



C b (7) 



r(i-7)x(7) 

r(l + 7)[-7X'(7)] 

f ,^ r j ^'(i)-v'(i--7) 



We now use the DIS factorization scheme, which is defined by the requirement for 
the Wilson coefficients Cf = 0, C\ = 1. Then the quark anomalous dimension simply 
becomes 

2iV /7 £ IS = h 2 ( lL )R™( lL ) (9) 

(cf. Eq. (D)). Fortunately, the expression for R® ls remains the same (Eq. (|[)), and as jl 
is per definition scheme independent, we can obtain fully resummed expressions in the 
gluon and quark sector. We get ^ qq from the following relation (jnl = Iqg)'- 

7 - = % { inl - l r 4) • (io) 

Explicit a s /n series expansions of the used terms can be found for example in |y| |l£ 



With these corrections jl, ^nl the singlet small x anomalous dimension matrix can 
then be written as follows: 

K^O + (™'*T^ f, T0 +o (<?)')■ (11) 

where 7,5, j v denote the yet unknown next-to- leading order gluonic contributions. Without 
them we will have a residual scheme dependence in the physical observables. Calculations 



in the Qo-scheme |23| suggest that the effects of setting them to zero will be comparatively 
small due to the dominance of 7 g9 and running coupling effects |24j]. Note that the non- 
singlet part of the renormalization group equation is not affected by small x contributions 
as singular as Jl, 'Jnl- 

3 Gluon based approach 

After this introduction to the formalism, we will focus on its implications on the descrip- 
tion of the measured structure function data. Our first example is the approach developed 



by Forshaw et al. in [16|]. The basic idea is to use only the LO small x contributions jl 



in the anomalous dimension matrix ( p]) and to neglect the quark singlet: 

/ r Q2 dk 2 

rin,Q 2 ) = f 9 (n,Ql)exp[ j ^ — lL \ (12) 



Following the argumentation of the previous section, we can include the renormaliza- 
tion factor R n into the gluon distribution before calculating structure functions: 

G(n, Q 2 ) = R n f 9 (n, Q 2 ) = R n p(n, Q 2 ) exp[Z n (Q 2 , Q 2 )) (13) 

with 

r Q2 rlk 2 

z n (Q 2 ,Ql)= %il. (14) 

JQI 



2 k 2 



If we want to transform this back into x-space, we have to perform the following 
complex integral: 

G(x,Q 2 ) = ^-j dnx- n R n f 9 (n,Ql)exp[Z n (Q 2 ,Q 2 )]. (15) 



If now the gluon input f 9 (x, Qq) is chosen as a simple step function JL6 



/"(re, Ql) = Af9(x -x)=* f 9 (n, Q 2 ) = Af^, (16) 

ft 

with a suitably chosen xq, there is no singularity other than n — 0, and the integration 
contour in Eq. ( |15| ) can be chosen as a circle around the origin. Then it is possible to find 
an analytic solution to the integral as a series of Bessel functions [|l( 



In order to form the structure function F 2 , one further integration according to Eqs. (|5]) 
and (^ has to be done. Forshaw et al. proposed to include the derivative of the coefficient 
function dC^/dlnQ 2 upon integration, although it is formally subleading in Eq. ([5]), since 
it may lead to contributions which are not negligible. The missing input 1*2(2, Qq) has 
been parametrized as A + Bx~ x , introducing three additional parameters. In |16|], xq 



was set to 0.1, and a (remarkably) low Aq^ d of 115 MeV was chosen. This leaves five 
parameters (Ql, Af, A, B, A) free to be fitted against structure function data. Using 1993 
data [0] below x = 0.01, they obtain a very good \ 2 i which is slightly better when the 
subleading terms mentioned above are included. It was our intention to test whether this 
approach still works for more recent data, which are more precise, and which especially 
extend down to lower Q 2 . 



Ql [GeV 2 ] 


data points 


M 


A 


B 


A 


x 2 


0.8 


257 


0.74 
0.41 


0.21 
0.42 


0.65 -10- 3 
-0.024 


-0.94 
-0.71 


650 
900 


without subl. terms 


1.0 


250 


0.80 
0.51 


0.23 
0.47 


0.0018 
-0.0066 


-0.79 
-0.57 


478 
646 


without subl. terms 


2.0 


232 


1.1 
0.71 


-0.032 
0.44 


0.12 
0.48 • 10~ 20 


-0.26 
-4.68 


337 
307 


without subl. terms 


2.0 


232 


0.39 


-2.40 


2.02 


-0.076 


1502 



(1) 



Table 1: Parameters for various fits to F 2 data. In the last row Aq^ d = 177 MeV 



3.1 Fits using the analytic formulae 

Our first observation is that the excellent agreement with the data in jl6| does not last 
very long. The analytic formula implies that in order to describe the low Q 2 data, one 
has to choose at least an equally low Qq. This already worsens the \ 2 f° r the 1993 data 
slightly (but it is still acceptable). If we now include the 1994 HERA data and the E665 
data p5|, the situation gets much worse, see table [l] for details. As in ||16| , we set xq = 0.1, 
and only used data with x < 0.01. One can see that the fits including the subleading 
terms are in principle better than those without. Nevertheless, the \ 2 °f & U these fits is 
not very convincing. It is better, the higher we chose Ql, so we focus in the following on 
Ql = 2 GeV 2 . 



As is visible in Fig. [I], the agreement with the data is not very good. The inclusion of 
the subleading terms (solid line) leads to curves which are "bent" down to the negative 
region for very small x. It is obvious that for these values of x this approach is not 



reliable anymore, as was already seen in ||16fl , too. With the 1993 data, this was not 
an issue, as there were no data in this kinematic region. But with the new data, this 
becomes a problem. Regarding the other curve without the subleading terms (dashed 
line), agreement with data is better, but we see also that the influence of F 2 (Q 2 ) ) is bigger: 



The gluon normalization M is smaller, and the parameters B and A acquire almost absurd 
values. 

The next question which arises concerns Aqcd- It is clear that this has a serious impact 
on the calculation, since the value of a s controls jl and thus the steepness of the gluon 
and the structure function. We tried to use a more realistic value of A.q CB = 177 MeV 
| 26| . This leads to an even stronger disagreement with the data. Moreover, relative to 
the other curves, the influence of the input i^x, Qq) is stronger, M being even smaller. 
This leads to the conclusion that the fitting procedure tries to reduce the influence of 
the (BFKL driven) gluon. To summarize, this simple analytic approach is disfavoured 
by the recent data. Note again that these data can be well described by a conventional 
renormalization group analysis. 

However, there remains one way out: If we abandon the simple step function like input 
gluon density, then the situation may be improved. It was noted in |fL6|l , that the analytic 
formula could be extended to more complex input densities. Unfortunately, we find this 
not to be the case: In the first place, one has to choose an ansatz whose Mellin transform 
has a finite number of n-plane singularities. This means e. g. for a standard gluon input 
of the form 

f(x,Ql)~x a (l-xf, (17) 

which transforms into a Beta-function, that either a or (3 have to be integer. While this 
is not too strong a constraint, the next problem arises immediately: For the analytic 
inversion it is essential that the input can be written as a power series in 1/n. However, 
the Laurent series stemming from Eq. (|T7| ) does not converge. Apart from that, one 
has to take into account that the expense of performing a fit based on a series of Bessel 
functions increases quickly, due to the nested sums involved and the greater number of 
parameters. So the analytic calculation is only feasible for the most simple form of the 
input gluon. We have chosen to simply resort to numerical methods, see below, but a 
more sophisticated analytic treatment is also possible p7 |. 



3.2 LO gluon fit 

Before any definite conclusions can be drawn on the applicability of this approach of 
Forshaw et al. 16|j , we have to overcome the restrictions mentioned above. We do this 



by transforming the moment space expression according to Eq. (|T5|) back to x-space 
numerically. This allows us to use a conventional form for the gluon which in x-space can 
be written as 

xg(x,Ql) = Nx- x (l + T 1 x)(l-xf. (18) 

The resummed gluon is evolved to Q 2 by multiplying this input distribution with R n 
and Z n (Q 2 ,Ql) according to Eq. (|T3|). An additional advantage of this method is that 
predictions for F 2 and Fl are obtained simply by multiplying in n-space the appropriate 
Wilson coefficients with the gluon solution before transforming to x-space numerically. 

We can now also exploit the fact that Z n can be calculated analytically. First note 
that due to the LO BFKL characteristic equation (^) we can write the measure of the 
integral as 

dlnQ =- lL 19 

Po n dji 

and thus we can rewrite the integral as 

r»lL(n,Q 2 ) ^ 

lL{n,Ql) 



z n (Q 2 ,Ql)= / d 11 - X {l) = lL 



n 

2 7E + — 



^■^ " 



T(1-7l) 



(20) 



This equation can be evaluated numerically once jl is known. 



We have tested a large set of fits using the gluon parametrization ( |18D shown above. 
The starting scale Qq was also fitted, but usually this parameter ended up at the upper 
limit set to the lowest Q 2 of the fitted points. We varied this lowest Q 2 from 2.0 GeV 2 
to 3.5 GeV 2 . Also we tested the effects of including the subleading term, of constraining 
the small x growth of the background F 2 ~ ^-o- 08 anc } f constraining the gluon power 
xg ~ x~ x to < A < 2. The results of the fits to the same data as in the previous section 
are summarized in Table 0. 



dCg 

dlnQ 2 

added 


rpbckgr. 
r 2 

~ X -0.08 


xg ~ x A 
0< A<2 


2.0 GeV 2 


X 2 /d.o.f. for Q 
2.5 GeV 2 2.8 GeV 2 


2 > 
3.0 GeV 2 


3.5 GeV 2 








3.08 


1.72 


1.35 


1.23 


1.00 






• 


3.12 


1.77 


1.41 


1.26 


1.04 




• 




3.17 


1.74 


1.43 


1.26 


1.03 




• 


• 


3.19 


1.78 


1.48 


1.29 


1.06 


• 






2.30 


1.02 


.679 


.625 


.583 


• 




• 


2.30 


1.02 


.680 


.626 


.583 


• 


• 




2.43 


1.33 


1.00 


.814 


.689 


• 


• 


• 


2.98 


1.33 


1.00 


.814 


.689 



Table 2: Fits using the formalism of Forshaw et al. and the general form of the gluon 



(|j~8|). The left part of the table shows the type of fit performed and the right part shows 
the x 2 P er degree of freedom for different lowest Q 2 of the fitted data. 



The fits including the subleading terms dC g /d In Q 2 describe the data much better. In 
spite of this we do not believe that these terms should be included, since they always lead 
to an obviously wrong behaviour in the small x region at low to medium Q 2 . A typical 
example is given by the dot-dashed curve in Fig. |^ that also displays other fits of Table 0. 
It shows the 2.8 GeV 2 fit of the bottom row in Table |2|. The negative contribution of 
these terms overwhelms the growth of F2 in the region ~ 3 — 10 GeV 2 at small x. The 
data constrain the fit, but it starts to fall right after the data point smallest in x. So 
the negative contribution on the one hand improves the description of the existing data 
due to a slowed growth, but on the other hand renders any prediction for future data at 
smaller x impossible. 



All the other curves in Fig. [| correspond to fits without subleading terms, but con- 
strained in F2 and xg (corresponding to row four in Table ||). For all fits except for the 
one with Q 2 > 3.5 GeV 2 , the fitted Ql equals the lowest Q 2 bin, so that the corresponding 
curve in that bin shows the F2 background. The fitted background of the exception is 
shown separately in the 3.5 GeV 2 bin. 



A general trend is that the fits get worse as lower Q 2 data are included. This is 
primarily due to the strong growth induced by the BFKL resummed anomalous dimension 
at low Q 2 . As is obvious by comparing the solid curves at 2.0 GeV 2 and at 2.8 GeV 2 
in Fig. |2|, the growth at small x is enormous over a short evolution length. Thus the fit 
forces the background to fall towards smaller x to be able to describe data at higher Q 2 
at all. To a lesser extent this is also true for the fits starting at 2.8 and 3.5 GeV 2 . The fit 
stays below the data at low Q 2 and small x in order to describe the bulk of data at higher 
Q 2 . Above 12 GeV 2 all fits describe the current small x data well. But if the precision 
and depth in x of the HERA data continues to improve, this limit may go up. 

We can conclude that despite some improvement upon using a more general gluon 
input, the fits still cannot describe the data well if started from low Q 2 scales. Even 
though the sub leading terms lower the \ 2 i a closer look reveals that their influence at 
small x leads to unphysical results. A word of caution has to be said about the influence 
of Aqcd- Lowering it to values as low as the one considered in the original paper [|1(J 
can halve the x 2 °f the fit, since the effects of the low starting scales are compensated 
somewhat. But in adopting such a low value of Aqcd one basically loses the connection 
to the large x region completely. If one attempts to vary Aqcd in fits from low starting 
scales Ql, then one sees soon that Aqcd always drops to the lowest limit set. This would 
of course not happen if there were large x data constraining the fit. So we here have 
chosen to fix it at a realistic value instead. 

4 Inclusion of the quark sector 

We have seen that all attempts to describe the structure function data using only a small 
x modified gluon density have failed. Thus the next logical step is to incorporate the 
more or less neglected quark densities. At the same time, we want to stay as close to the 
successful conventional renormalization group equations as possible. A method to achieve 



this was suggested by Catani et al. |14] and used in the paper of Ellis et al. [|17]], but only 



10 



by evolving an existing set of parton distributions (MRSDq, p9|). In the following, we 



describe this approach briefly and present the results of our fits based on this method. 

The basic idea is to take the renormalization group equations in their two-loop-form, 
and to modify the anomalous dimension matrix according to Eq. (|TTD . The importance of 
implementing the fundamental energy-momentum conservation for the modified equations 



has already been stressed in [17]. Energy- momentum is of course always conserved in 
the conventional formalism. In Mellin space this means that the first moments of the 
anomalous dimensions have to vanish: 

£ 7y (n = l) = 0. (21) 

i 

The all order small x resummations violate this equation and it has to be enforced by 
hand somehow. Additionally, when we combine the two-loop with the small x expressions, 
we have to avoid double counting of the leading terms which appear in both the two-loop 
expressions and the resummed expressions. This can easily be done by simply subtracting 
the first term in jl, and the first two terms in Jnl- As we are now working in next-to- 
leading order, we have to choose a definite factorization scheme. In the following we will 
work in the DIS factorization scheme, since in this scheme there exists a fully resummed 
expression for jnl, as mentioned earlier. 

Concerning energy-momentum conservation, several ansaetze have been suggested in 
the literature. Ellis et al. Jl7j used two different ones, a "hard" one where the small x cor- 
rections 7jj(n) are replaced by Jij(n) —7^(1), and a "soft" one where these corrections are 
multiplied with a factor (1 — n). Blumlein et al. [|TBJ additionally proposed the "conserving 
factors" (1 — n) 2 and (1 — 2n + n 3 ). Although these are all arbritrary implementations of 
the momentum sum rule, one can study with these different factors the possible impact 
of yet unknown higher order corrections^. In Jl8| it was already shown that with these 



factors it is possible to suppress the small x corrections, or even to overcompensate the 



1 Ball and Forte |19| ] have circumvented the energy- momentum conservation problem in a different 
way. They used a factorization scheme in which the unknown higher order terms 75, 7^ in the anomalous 
dimension matrix are defined by the momentum sum rule Eq. (plf). 

11 



expected growth at small x. 

There are some theoretical problems associated with these ansaetze. First, we consider 
the hard implementation. While this does not modify the n-dependence of the resummed 
corrections at all, it requires their evaluation at n — 1. This causes a problem if the 
starting scale of the evolution is low, since the resummed n-plane pole n r = o7 s 41n2 can 
in this case pass through n = 1, making further calculations impossible. For example, the 
MRS Rl fit U uses A QC d = 241 MeV and a starting scale Ql of 1 GeV 2 , leading to n r = 1 
at Q 2 = 1.07 GeV 2 using the two-loop formula for a s . Moreover, the generated structure 
functions are still much too step to describe the data, so we will not use it further. 

Second, factors involving higher powers of n than n 2 (like one of those in [18|]) cause 
a factorization problem. If we remember that the leading term of 7 99 beyond the two- 
loop expression is proportional to (a s /n) 2 , such conserving factors imply that j^l times 
the conserving factor does not vanish anymore for n — > oo. This means it cannot be 
inverted to x-space on its own. Since the anomalous dimensions are folded with parton 
distributions which for reasonable choices of the input fall off better than 1/n for n — > oo, 
this is not a problem for the evolution program, but it is still an indication that this 
product is not correctly factorized anymore. For this reason we do not examine even 
higher powers of n in the conserving factor. 

Now we have to comment on our implementation of the resummed expressions into the 
two-loop evolution program. Since we are using a program which is based on the two-loop 
analytical solution in moment space of the renormalization group equation || |3(J, this 
is not as simple as it would be e. g. in a Runge-Kutta based x-space evolution program. 
Let us first write down the singlet renormalization group equation in the two-loop form 
including the small x resummations: 

df(n,Q 2 ) 



T (0) (n)+(g) 7 (1) W+7 re >,«.) 3-x^/M 2 ) (22) 



a s . fn\ , x fa 

27 



~4w a s (47r) 2Ct s 



da s 
The vector / consists of the quark singlet f s and the gluon density f 9 , and we employed 
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the two- loop approximate solution for a s : 

a s (Q 2 ) 1 



p\hiln(Q 2 /A 2 

17" ~/? ln(Q7A 2 ) /3 3 [ln(Q 2 /A 2 )p 



(23) 



The solution of Eq. ( p2|) is not trivial, since we are dealing with a matrix equation; 
furthermore, it is unclear how to incorporate the resummed expressions involving higher 
powers of a s into an approximate solution which is accurate only to 0(a 2 ). We do this 
by separating the treatment of higher order terms ot\, i > 3 from the resummed terms 
a s (a s /n) k , k > 2: 



f(n,Q 2 ) 



exp 



-7(°)(n)ln^3 
/5o 7 W a s (Ql) 






(24) 



as(Q 2 ) 



da 



f es (n,a s 



o.Wg) 



*q,2 i ft Q 3 



+ o q,o: 



a. 



n 



with k > 2. 

Since 7 res starts with a s (a s /n) 2 after subtraction to avoid double counting, the leading 
power of the integral is (a s /n) 2 . In expanding the matrix exponential we would encounter 
commutators of all three terms. But the commutator of the conventional NLO term with 
the resummation integral is suppressed by a s . This contribution would be of the same 
order as the commutator of the error of the resummation with the LO RG term and can 
be ignored. Thus we can attach all the resummed terms to 7^ via 



7 (1,rcs W s 



7 (1) (™) + 



a s (Q 2 ) - a s (Q 2 ) 



da ■ 



7 res (ra, a) 



@Q_rv 2 -I- & 1 <r 3 ' 



(25) 



For the final solution we then treat the resummed anomalous dimension 'y' 1 ' res ) exactly 
like the original 7^ 1 - ) . The exponential of the matrices is expanded and terms of 0(a 2 ) not 
stemming from the resummation are dropped. Thus, concerning the resummed terms, 
the solution is accurate up to 0(a 2 (a s /n) k ) consistent with the error inherent to the 
resummation itself. 
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We have checked that using this procedure we can reproduce the results of |H]]. But 
we want to improve on the results shown there by on the one hand refitting the input 
parton distributions and on the other hand by lowering the starting scale Qq to values 
used in recent conventional RGE fits. We expect that lowering the starting scale leads to 
the same kind of problems as encountered in the previous section, namely a steep growth 
in the region of small x. Thus it is necessary to refit the parton distributions to test if 
one can compensate this growth. 

In principle one should do a fit to all relevant data to constrain the (presumably) 
universal new input parton densities and a s as much as possible. Then the analysis 
would be truly competitive to conventional RGE analyses like || [|. But the small x 
resummations for most of the processes are unknown. Furthermore the computing time 
would be prohibitive due to the more complicated expressions. Even if we just concentrate 
on F 2 , the calculations are already rather time consuming. 

Thus we take the following approach: we use the optimal MRS Rl conventional RGE 
fit || as our starting point. The ansaetze for the input parton distributions can be found 
there. We assume that at large x the new parton distributions would be similar to the Rl 
ones, since the resummation effects should be small at large x and since there are many 
experiments constraining the parton distributions in this region. In practice we calculate 
the original Rl partons and F 2 at x = 0.05, 0.15, 0.4 and compare them to the new fit 
using an artificial error of one percent. This assumed small error is necessary to force the 
fit to match the large x region, since as we will see the small x region is not well described 
and would overwhelm any weak constraint at large x. The valence quark distributions are 
not directly affected by the small x resummations. They could be influenced indirectly 
by trading momentum with the quark sea, but they are well constrained by experiments 
and negligible in the region of small x. In order to save computing time we do not fit the 
parameters of the valence quark input, but keep them fixed at the Rl values. We also 
keep the value of Aqcd fixed, since it is mainly constrained by experiments in the large x 
region as discussed in ||. 
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1UDCD 

We allow for an additional degree of freedom by fitting Aa = —0.3 of xA = x(d — 
u) ~ x~ Xa . The correct u — d flavour symmetry breaking is taken care of by the large x 
constraints of the sea. The glue and the sea are now fitted at x < 5 • 10~ 2 to the data. 
Of course in this region the HERA data |l], |2|, |3| are dominant. Fitting to the whole bulk 
of HERA data would again be too time consuming. So instead we use the Rl fit, which 
describes HERA data very well, to generate F 2 data points in the x range that is covered 
by HERA. The errors are adjusted to match those of the experiment. Finally it should 
be mentioned that the appropriate adjustments of the parton distributions due to the 
scheme change MS — > DIS rcs have to be madeQ. This is done automatically at the starting 
scale Qq, so that only the input distributions are affected. Since we are mainly interested 
in how the parameters are changed when the small x resummations are switched on and 
since we keep the valence quarks fixed, only the conventional Wilson coefficients are used 
for the transformation. This does not limit the fits, but the obtained parameters should 
not be used for a direct MS calculation. 



Let us mention that in the previous fits by Ball and Forte |15, |2(J only the small 
x powers ~ x~ x of the glue and sea and the starting scale Ql were treated as free pa- 
rameters. To be precise, they evolved several sets of existing parton distributions, which 
are compatible with large x data, to Ql using conventional two loop RGE. The original 
ansatz was then refitted to these resulting partons. Finally a fit to the experimental data 
using these new input parton densities was performed, in which only the mentioned three 
parameters were allowed to vary. The large x behaviour of the partons can be reproduced 
by rearranging the various large x parameters, so that for example the normalization can 
change. This does of course affect the small x region as well, so that we instead find it 
necessary to fit the complete ansatz. Also they introduce a "reference value" xq above 
which all resummed expressions are switched off. Although their motivation to separate 
the large from the small x region is clear, we see no need in the formalism to introduce 
such an arbitrary scale. Nevertheless, their conclusions are quite similar to ours, since 



2 The superscript "res" merely signifies the inclusion of the small x resummations, not a change of the 
scheme definition. 
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in a later paper [^UJ xo is found to be almost zero, thus excluding any need for small x 
resummed expressions in the evolution equations. So we will basically confirm their re- 
sults with our different method, which determines all parton parameters for the small x 
modified evolution by directly fitting to experimental data. 

To see the influence of the resummation on the parton densities we show in Fig. ^] 
the ratio of the resummed gluon and quark singlet to the standard one. As expected the 
ratio goes to one for large x, but the expected slower convergence of the sea curves is also 
clearly visible. For larger Q 2 the influence of the resummation is much reduced. Most 
of the difference is picked up at low Q 2 values which is understandable considering the 
moving resummed pole. To test the impact of yet unknown higher order contributions 
as mimicked by our conserving factors, we are especially interested in the first term ~ n; 
but changing its coefficient requires the inclusion of higher powers of n. Comparing the 
(1 — 2n + n 2 ) and (1 — 2n + n 3 ) curves, we see that their influence is strong. In the latter 
case, the partons are even smaller than the standard ones and a fit using this factor leads 
to an unacceptable fall of F-z in the small x region. Thus it is evident that present small 
x resummations have no predictive power. 

On the other hand we have already mentioned that exactly the term ~ n 3 and all 
even higher powers spoil the transformation back to x space. If we require the anomalous 
dimensions to be transformable, the maximum power allowed for the conserving factor is 
n 2 . In combination with the strong large x constraints we then always get partons and F 2 
larger than the ones obtained by conventional NLO evolution. We checked this by fitting 
the conserving factor as well, i. e. we fitted a of (1 — a ■ n + b ■ n 2 ) with b = a — 1. The 
fit gives a = 1.2944 and satisfies the large x constraints well, but the (1 — n) 2 curve does 
better in the small x region. Actually if we only fit a and leave all other parameters at 
their original MRS Rl value, we get a = 1.2651. So for these values of a the effect of the 
small x resummations at large x seems to be minimized. 

Such large effects from subleading terms are caused by the resummed BFKL pole. For 
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Qq = 1 GeV 2 it lies around n r ~ 1. So considering a contour close to the resummed 
pole we see that the main contribution of the resummed pole comes for regions of n 
where the change n 2 — > n 3 in the conserving factor means a considerable change of the 
subleading terms. Of course this depends on the starting scale Ql as well, for example at 
Qo = 4 GeV 2 we have n r ~ 0.7. So on top of the change of the resummation itself the 
relative importance of the subleading terms will change. This arbitrariness is just a sign 
of our ignorance concerning the subleading terms, but at least we can exclude powers of n 
higher than two from the conserving factor to avoid picking up terms that will not occur 
in the real subleading expressions. 

The results for the partons translate directly to F2(x, Q 2 ) shown in Fig. |j. As expected 
the resummed calculations always overshoot the data at small x. The influence of the 
terms down by one power in n is strong as can be seen by the spread in the curves. The 
"dips" in the curves are somewhat artificial, basically they are due to the too strong rise 
of the calculated F 2 compared to the data. Then y 2 is minimized by adjusting the partons 
so that the curves are below the data at larger x first. Even the curve using the (1 — n) 2 
factor is too steep at small x; this is even more true for the other curves. The curve with 
the fitted conserving factor a = 1.2944 mainly leads to an optimal large x behaviour. It is 
interesting to note that if we leave out all large x constraints and fit the conserving factor, 
we get a good description (x 2 /d.o.f . ~ 1) of the small x data for a = 2. But then at larger 
x we are totally off the F 2 and parton constraints. Due to the coincidence that the fitted 
a for the best small x description is close to the (1 — n) 2 we tested for the complete fit, 
we can claim that the (1 — n) 2 curve is basically the best if one wants to describe all F 2 
data but emphasizes the small x region. 

A problem of all fits has to be mentioned. If we leave Aqcd as a free parameter it 
always drops to low values. Despite the still large errors on Aqcd and despite the fact 
that in principle all determinations of a s should incorporate small x resummations to be 
consistent, we believe that Aq^ ~ 200 MeV should be taken as typical for a "low" a s . 



Experimentally low values of a s stem from data at large x [31] where resummation effects 
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x[g\S\A](x 1 Ql) = 


RGE 


Resummed 


Ax~\l - x) v (l + e^fx + 72;) 


MRS Rl 


{l-nf 


a = 1.2944 


(A) 


24.5 


0.94777 


36.09 


\ 


-0.41 


0.26536 


-0.547 


Gluon rjg 


6.54 


3.9801 


6.866 


e 9 


-4.64 


-1.5297 


-4.483 


7 9 


6.55 


4.2279 


6.135 


A s 


0.42 


0.41647 


4.801 


As 


0.14 


0.036837 


-0.3928 


Vs 


9.04 


10.873 


8.815 


Sea €3 


1.11 


5.1104 


-0.330 


(e A = 0, i] A = r]s) 7s 


15.5 


22.073 


0.441 


A A 


0.39 


0.084878 


0.001299 


A A 


-0.30 


-0.59879 


-0.1086 


7a 


64.9 


72.461 


1473 


^-(F 2 ),x<5-10- 2 




1.96 


12.0 


^y(F 2 , partons), x > 5 • 10~ 2 




126 


6.04 



Table 3: The fitted partons using RGE plus small x resummations in comparison with 
the original RGE MRS Rl partons [§. 

should be small, so we expect that this holds true in the cases considered here. The effect 
of lowering Aqcd is of course greatest for the largest resummation effects, so we show a 
curve with the conserving factor (1 — n) and Aq^ = 200 MeV. The x 2 /d.o.f. at small 
x is for the lower A QCD reduced from 21 to 11. A straight fit would yield Aq^, =180 
MeV, but this value should not be taken too seriously, since it depends on how strictly 
we implement the large x constraints. 



In Table |3] we show the parameters of those two fits that can be considered optimal. 

n f =4 



A g is calculated via the momentum sum rule and Aq CD 



241 MeV. S = 2(u + d + s) 
is the total sea quark distribution, and A = d — u; for further details and the valence 
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distributions refer to ||. The x 2 for the (1 — n) 2 conserving factor is dominated by the 
partons constrained by an artificial one percent error. A more realistic choice of a five to 
ten times bigger error would of course reduce the x 2 considerably. But we see that even at 
small x the quality of the conventional RGE fit is not matched. Most of the x 2 comes here 
from the points with smallest x at the lower Q 2 bins. Looking at the parameters we see 
that the input sea has turned flat whereas the input gluon is now growing with smaller x 
instead of being valence like. The normalization of the gluon is lowered correspondingly. 

The input sea of the a = 1.2944 fit has become valence like. The gluon has stayed 
more or less the same. On the one hand we see that the fit cannot describe the small x 
data at all. On the other hand at large x the MRS Rl partons and Fi are reproduced 
well. We can be sure that this fit would be compatible with experimental data at large 
x. Due to the necessary artificial constraint at large x it makes no sense to add the x 2 
for all x, so one cannot conclude that the a = 1.2944 is the overall better fit. A more 
sophisticated treatment of the large x region would be needed to determine a true best 
fit. But since there are still parts of the NLO small x resummation missing and since the 
conserving factors give just an estimate of the influence of subleading terms, we feel that 
not much could be gained by improving the large x treatment for the time being. 

In order to minimize the systematic error introduced by the missing terms, we also 



employed the Q -scheme p3[. As already mentioned, at least the qq contribution to the 



missing gluonic small x anomalous dimensions is small |24] in this scheme. The Wilson 



coefficients for F 2 and the small x anomalous dimensions stay the same as in the DIS 
scheme, except for setting R n — > 1 in Eq. (^|). Since R n = 1 + 0[(a s /n) 3 }, no changes are 
introduced in the conventional two loop expressions. The effect on F2 of setting R n — > 1 
is explored in Fig. ||. There the deviation (F 2 — -F 2 DIS )/F 2 DIS is displayed for different Q 2 
starting from the same input parton distributions. The conserving factor (1 — n) 2 and 
the parameters of the corresponding fit of Table |3| were used for the evolution in both 
schemes. 
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The deviation is less than ten percent in the x and Q 2 range of the data considered 
here. Also we can see from Fig. [5] that the naive expectation that the DIS sea and thus F 2 
should be larger at small x due to R n > 1 only becomes true for x values at the edge of 
current data. We can even expect a Qo-scheme fit to fare worse in the x < 5 • 10~ 2 range, 
since the bulk of HERA data is in the region where the Qo-scheme calculation results in 
an even larger F 2 . We have confirmed this expectation by repeating the (1 — n) 2 fit in the 
Qo-scheme. Actually the growth of F 2 in the main data range is enhanced by the strict 
large x constraints. This is shown by the dotted line in Fig. ^ which display the deviation 
at 65 GeV for the partons fitted in the Qo-scheme. The strong large x constraints force 
the deviation to zero at three points, which leads to the notable, artificial oscillation in 
that region. Thus the resulting large \ 2 a t small x should only be taken as indication that 
switching to the Qo-scheme does not give us a better estimate of the missing contributions 
at the moment. 

5 General Discussion and Conclusions 

Figure ||| presents a different look at our results. Here we show A = d\nF 2 /d\n(l/x) as 
determined by the HI experiment || in comparison with theory. To obtain A a fit of 
the type C x~ x is made for x < 0.1 for the HI data. This should be understood as a 
cut, so that the actual highest and lowest x value used depends on the data. The two 
conventional NLO RGE calculations MRS Rl |§ and GRV '94 [§2j] have been treated by 
us in a similar way using the same highest x values but extending the lowest x to 10 -5 . 
We see the perfect match of the MRS Rl fit, and show for comparison the dynamical 
GRV '94 [[32] predictions which, however, are sensitive to the precise choice of the input 



scale. 

We also show in Fig. [] the slope predicted by the pure LO BFKL evolution. For 
this we simply used the A of our parametrization [ lQf . For the "resummed" calculations 



presented in the last sections we again used the fitting method to obtain A. But to avoid 
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the "dips" of the curves presented in the last section we use a x < 1CT 3 cut and for the 
ones of the section |3] we even have to use x < 3 ■ 10 -4 . We show in Fig. |6] A for all starting 
scales Qq of the fits displayed in Fig. |2| and for the different fits displayed in Fig. |j. 

We see that the big spread in the predicted A for all the BFKL inspired methods is 
diminishing towards higher Q 2 . But it lies always above the data except perhaps at very 
high Q 2 . So the resummed F 2 is always growing too fast in the small x region, since 
at very large Q 2 the conventional RGE terms are expected to become dominant. This 
is also probably the explanation why the pure LO BFKL A is the only one that keeps 
growing. There are no conventional RGE terms that can take over at high Q 2 for this 
curve. The fact that the LO resummed calculations give predictions for A that are in the 
same ballpark as those including additionally the NLO resummations can be attributed 
to the dominance of the BFKL pole. 



The failure of the pure LO BFKL formalism |10|, |TT| to describe the experimental F 2 



data is not so surprising. There has never been any guarantee that the resummation of 
only the leading logarithms in x would be appropriate in the kinematic region explored 
by current experiments. The inclusion of a larger part of the total contributions of per- 
turbative QCD should naively lead to a better description of the data. Our calculations 
using the methods of Forshaw et al. show us clearly that it is not sufficient to use just 
gluons and the LO resummed anomalous dimensions in RGE type calculations. The most 
promising method should be to use all that is known: the NLO RGE and the LO and part 
of the NLO resummations in x. In doing so we immediately encounter the first hint of 
trouble. The resummations in x always violate the fundamental energy-momentum sum 
rule. 

Thus we can really only rely on the calculations if the results are not changed much 
when we implement this constraint in various ways. In fact the results depend strongly 
on the implementation and one cannot give reliable predictions for F 2 at small x at all. 
The main effect stems from rather large n and thus variations in the conserving factor 
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lead to changes that are formally subleading but numerically important. If we limit our 
calculations to conserving factors that do not introduce terms absent from higher order 
anomalous dimensions, then we can still see a trend in the predictions: again F% grows 
too strongly at small x. 

Our results could be taken as an indication that the still missing NLO resummation 
pieces will not be sufficient to improve the stability of the calculation and its outcome. 
This would prove that the small x resummations do not lead to a stable perturbative 
series in contrast to the RG calculations. Of course this is entirely possible, but it would 
question all future work on this subject. But there are some important caveats concerning 
such a strong conclusion. 

A strong growth that can not be suppressed by adjusting the parton distributions is 
only encountered if the starting scale Ql is below approximately 3 GeV . To a lesser 
extent the variations induced by different conserving factors increase for lower starting 
scales. Also the fits can always be improved significantly by lowering Aqcd- Only a true 
fit to all the large x data would tell us if A QCD can be determined uniquely and if the 
small x region is well described with the value obtained. But our fits already indicate 
that probably this would not be the case. These problems are all connected to the strong 
BFKL pole, since the fits always improve if a s becomes smaller and the pole moves to the 
left in the n-plane. 

It has to be mentioned as well that some closely related calculations do not seem to 
encounter these problems. The "physical anomalous dimensions" calculations |H], |22| re- 
late only physical observables F 2 and F^ and thus any factorization scheme dependence is 
avoided. This fit to the data even seems to be preferred over conventional RG calculations. 
We note that due to the relative order Fl ~ a s F2 such calculations should be limited to 
leading logarithms until the next-to-next-to-leading piece of the coefficient functions of 
Fl is determined. This fit also implies a rather low (LO, 4 flavours) Aqcd= 100 MeV. 
The input scale of the evolution Qq is determined as 40 GeV 2 and the fit becomes un- 
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competitive for Q 2 , > 3 GeV 2 . Finally a new scale All is introduced at which the "inputs 
become nonperturbative" |[21|| . The usual choice All = Qo would simplify the formulae 
used considerably. But this would spoil the agreement with data, since Ql would be low. 

We suspect that if this kind of fit was redone starting from a low scale Qq = All, then 
the same kind of problems as we have encountered would occur. Calculations using the 
colour dipole model are also giving good agreement with recent F 2 HERA data [J3~3"| , but 
note that they only used a fixed strong coupling. Another successful approach J34j based 



on the CCFM equation |15], |35| , does not include the quark sector, uses approximations 
only valid at small x and represents a LO resummation. But the imposition of a purely 
kinematic constraint on the gluon ladder introduced in p4j suggests that higher order 
terms could suppress the growth in the gluon sector and thus effectively shift the BFKL 
pole to the left. 

For a final judgment on the small x resummations we will have to wait until the 
calculation of the NLO pieces will be completed. Also it is probable that with a starting 
scale of around 4 GeV 2 it will not be possible to rule out such contributions. Only a 
true fit to all the usual data using small x resummations consistently will tell us if the 
assumption that Aqcd can not be lowered too much because of the large x region is 
correct. Keeping this in mind, we suggest that currently available methods and data 
strongly favour conventional RG calculations. If this statement survives the test of time, 
it will mean that the small x resummations do not represent a good perturbative series. 
This would shed doubt on all calculations involving them. 
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Figure Captions 

Fig. [I] Structure function obtained with the analytic prescription of Forshaw et al. fl(| . 
with refitted parameters according to the shown data from HERA [fj, [3| and E665 
[25J]. The solid curve is calculated without the subleading term dC 2 /dln Q 2 , whereas 
the others curves include them, with two different choices of Aqcd- 

Fig. |2| Fits to F 2 data as in Fig. Q but using the gluon in Eq. (|T7D. The lower Q 2 cut 
of the fitted data is varied and the fitted starting scale Qq ends up at this value 
except for the fit Q 2 > 3.5GeV 2 . For this fit the background at Q^ = 3.36 GeV 2 is 
displayed in the 3.5 GeV 2 bin. The dot-dashed curve is an example of a fit including 
subleading terms. 

Fig. |3] Impact of the resummed terms in the evolution equation on the singlet (left) and 
gluon part (right) relative to conventional two-loop calculations. Also the influ- 
ence of different methods to restore energy-momentum conservation is shown. The 
unmodified MRS Rl distributions || are used as input. 

Fig. |I] The results for F 2 of fits using different conserving factors in a NLO calculation 
including small x resummations. The data used for the fit is shown as stars and 
experimental data as in Fig. p] is shown for comparison. 



Fig. |5| The deviation of F 2 calculated in the Qo-scheme [23|, |24| from the DIS scheme at 



several Q 2 . The parton distributions of Table |3| obtained in the (1 — n) 2 DIS scheme 
fit are used as input. The dotted line shows the deviation at 65 GeV 2 , when partons 
fitted in the Q - sc heme are used for F^ ■ 

Fig. |6] The slope of F 2 of the calculations using small x resummations compared with 
experimental data and two conventional NLO RGE calculations. Results of the 
same method are shown in one line style. The main parameter varied to obtain the 
different curves of one method is displayed at the curve. 
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